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Abstract 

The  local  stability  of  several  octahedral  species  A1F63-,  GaF63-, 
InF6^',  T1F63*,  TaF6‘,  ZrFe^',  LaFg^*,  and  LaCla^*  was  studied  using  ab 
initio  methods.  All  systems  were  found  to  have  minima  at  octahedral 
geometries.  The  electronic  stability  of  these  species  was  considered 
using  electron  propagator  theory  (EPT)  and  the  M0ller-Plesset  (MP) 
perturbation  scheme.  The  results  indicate  that  triply-charged 
molecular  species  may  be  electronically  stable  in  the  gas  phase.  Due 
to  the  slow  convergence  of  the  EPT  and  MP  series  and  limitation  in 
our  one-electron  basis  set,  the  final  conclusion  requires  a  more 
advanced  theoretical  treatment.  The  doubly-  (ZrF62-)  and  singly- 
charged  (TaFg-)  species  were  found  to  possess  very  high  electronic 
detachment  energies  of  5.0  and  10.6  eV,  respectively.  TaFe’  is 
predicted  to  be  thermodynamically  stable  and  ZrF^^-  slightly 
unstable  with  respect  to  the  unimolecular  decomposition  ML6"*  — > 
ML5(n-l)-  +  L-. 
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I.  Introduction 


Multiply-charged  anions  are  challenging  species  from  the  point 
of  view  of  electronic  structure  theory.  Coulomb  repulsion  between 
excess  negative  charges  has  to  be  counterbalanced  by  an 
advantageous  distribution  of  the  total  charge  to  make  the  system 
electronically  and  geometrically  stable.  The  multiply-charged 
systems  are  usually  thermodynamically  unstable  but  if  their 
lifetimes  are  long  enough,  they  may  be  considered  as  promising 
candidates  for  high  energy  density  materials.^ 

The  discovery  of  the  dianion  of  benzol  [cd]  pyrene-6-one  by 
Dougherty^  opened  the  field  of  doubly -charged  molecular  species  in 
the  gas  phase.  An  excellent  review  of  the  subject  from  both  the 
experimental  and  theoretical  perspective  has  recently  been 
published  by  Compton. ^  Briefly,  Leiter  et  al.  found  doubly  charged 
anions  of  the  oxygen  clusters  (02)ii^‘-‘^  Schauer  et  al.  reported  on  the 
long-lived  Cn^*  (n=7  to  28)  species,^  and  dianions  of  fullerenes  were 
found  as  well.®  Korobov  et  al.  estimated  the  electron  affinity  of  PtFe* 
to  be  positive^  and  Maas  and  Nibbering  have  produced  doubly 
charged  dicarboxylate  negative  ions.^ 

The  challenge  to  supersede  the  Coulomb  repulsion  between  two 
extra  electrons  provoked  theoretical  studies  as  well.  Gutsev  and 
Boldyrev  estimated  maximal  value  of  the  second  EA  of  hexafluorides 
MF6  in  the  region  of  1-2  eV  on  the  basis  of  the  electronic  properties 
of  the  F6  frame.^a  They  also  calculated  second  EA  (SEA)  of  PtFg  by 
DV-Xa  method  and  found  SEA=0.7  eV.^a  Miyoshi  and  Sakai  predicted 
positive  electron  affinities  for  CrFe*  and  MoFe'.^O  Recently,  Ewig  and 
van  Waser  found  CFe^-  to  be  geometrically  metastable  at  Oh 
symmetry.il  Smaller  tetrahedral  (BeF47*)  and  triangular  (LiFs^-, 
NaFs^-,  KF32-,  LiCla^',  NaCla^-  and  KCI32-)  dianions  consisting  of 
halogen-type  ligands  and  a  metal  atom  have  recently  been  carefully 
studied  by  Cederbaum  and  collaborators.  12-1^  They  presented 
convincing  evidence  for  the  stability  of  these  species  with  respect  to 
fragmentation  and  electron  ejection  as  well  as  a  kinetic  stability  of 
these  anions.  Although  these  dianions  are  not  stable 
thermochemically  with  respect  to  dissociation  into  two  singly 
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charged  anions,  all  these  species  have  a  high  barrier  and  therefore 
have  a  long  lifetime.  The  ionic  nature  of  these  systems  was 
recognized  as  a  key  factor  in  compensating  Coulomb  repulsion 
between  two  extra  electrons. *3. 14  Boldyrev  and  Simons  suggested 
TeFg^’  as  the  MXn^*  species  with  the  largest  electron  detachment 
energy,  15  because  delocalization  of  two  extra  electrons  through  eight 
electronegative  fluorines  occurs.  The  same  authors  also  searched  for 
linear  doubly-charged  anions  and  suggested  Mg2S32'  as  the  smallest 
linear  or  quasilinear  electronically  stable  species.  1^  Second  electron 
affinities  of  large  linearis  and  of  trigonal  planarize  carbon  cluster 
have  been  found  theoretically.  However,  some  free  doubly  charged 
anions  such  as  S04^*,  Se042'  and  Te042-*18  have  been 

predicted  to  be  electronically  unstable  nonwithstanding  that  these 
anions  are  postulated  in  solid  state,  in  solution  and  in  molten 
substances. 

As  far  as  the  existence  of  more  than  doubly  charged  anions  is 
concerned,  neither  theory  nor  experiment  has  delivered  a  definite 
answer  yet.  Early  reports  on  were  interpreted  as  harmonic 

artifacts.^  Miyoshi  and  Sakai  concluded  from  their  ab  initio  Cl 
calculations  that  AuF6  is  electronically  unstable  at  octahedral 
geometries lOb  and  Boldyrev  and  Simons  have  shown  that  POa^*, 
AsOa^*  and  SbOa^*  are  not  electronically  stable.!*  A  scheme  for 
building  multiply  charged  complex  anions  has  been  recently 
proposed  by  Scheller  and  Cederbaum^h,  and  a  similar  approach  was 
also  briefly  discussed  by  Gutsev  and  Boldyrev. 9 a 

In  this  work  we  explore  local  stability  of  singly-,  doubly-,  and 
triply -charged  closed-shell  octahedral  anions  MLe.  By  local  stability 
we  understand  existence  of  an  Oh -symmetry  geometrical  minimum 
at  which  the  value  of  the  electron  vertical  detachment  energy  (VDE) 
is  positive.  Of  course,  these  two  features  are  necessary  but  not 
sufficient  for  long-lived  (meta)stability.  The  latter  requires  that  the 
closed-shell  anion  remain  electronically  stable  for  a  wide  range  of 
molecular  geometries  which  are  probed  in  zero-point  vibrational 
motion  and  that  the  barrier  height  for  unimolecular  decomposition  of 
these  species  be  sufficiently  high  and  wide.  These  topics  are  beyond 
the  scope  of  the  present  paper.  They  will  be  addressed  in  the  future 


for  the  most  promising  candidates  selected  on  the  basis  of  their  local 
stability  at  octahedral  geometries. 

II.  Systems  of  Interest 

Octahedral  hexafluorides  represent  a  good  choice  for 
electronically  stable  multiply  charged  anions,  because  extra  charges 
are  distributed  through  the  large  number  of  electronegative  ligands 
and  repulsive  ligand-ligand  interactions  are  not  very  high.^O  Fully 
delocalization  of  extra  electrons  is  expected  for  MFk+r,  MFk+2^', 
MFk+a^'  and  for  MFk+n"*  in  general  cases,  where  k  is  the  maximal 
formal  valence  of  the  central  atom  M  and  n  is  the  negative  charge  of 
the  multiply  charged  anion.  The  central  atom  M  should  be  from  the 
group  V  of  the  Periodic  Table  (P-Bi,  V-Ta)  for  singly  charged  anions 
(k+l=6),  from  group  IV  (C-Pb,  Ti-Hf)  for  doubly  charged  anions,  and 
from  group  III  (B-Tl,  Sc -La)  for  triply  charged  anions.  For  ail  such 
systems,  octahedral  closed  shell  structures  are  expected. 

Quantum  chemical  calculations  (DVM-Xa)  of  the  nonmetallic 
hexafluorides  MF6'  and  hexachlorides  MCl6'  have  shown  high  vertical 
electron  detachment  energies  (VEDE)  for  these  species:  PFe*  (6.8  eV), 
AsF6*  (7.2  eV).  SbF6*  (7.0  eV),  PCle*  (4.3  eV),  AsCle*  (5.4  eV)  and 
SbCl6*  (5.6  eV).2l  While  experimental  electron  detachment  energies 
for  these  anions  are  not  known,  the  theoretical  results  agree  well 
with  the  thermochemical  stabilities  of  noble  gas  compounds  such  as 
KrF+[SbF6]-.22  Xe2F3+(AsF6]-’23  and  others  (see  Refs  24  and  25). 

The  doubly  negative  nonmetallic  hexaflouride  anion  CF62-  is 
geometrically  stable.  However,  the  electronic  stability  has  been 
studied  only  for  SiF62-,26  what  was  found  to  be  electronically  stable, 
with  a  vertical  electron  detachment  energy  equal  to  0.2  eV.  No 
experimental  data  is  available  for  the  stability  of  free,  doubly 
charged  MF62-  anions,  where  M=C-Sn. 

The  electronic  stability  of  triply  charged  MFk+3^*  anions  with 
nonmetallic  central  atoms  has  not  yet  been  reported. 

Unusually  high  first  electron  affinities  of  hexafluoride 
molecules  containing  transition  metal  atoms  were  recognized  long 
ago. 27.28, 29  For  instance,  the  electron  affinity  of  AuFe  is  believed  to 


be  ca.  10  eV.27.29  Theoretical  calculations  of  Gutsev  and  Boldyrev2  0 
confirm  a  high  electron  affinity  of  all  3d,  4d  and  5d  metal 
hexafluoride  species  but  also  provided  a  theoretical  explanation  for 
the  unusual  electron  affinities  of  VF6,  NbF6,  and  TaF6  as  well  as  for 
hexafluorides  of  the  metals  from  the  end  of  the  transition  period 
(CuF6.  AgFg,  AuF6,  see  Ref.  20  for  details).  Miyoshi  and  Sakai 
predicted  the  dianions  CrF6^*  and  MoF62*  to  be  electronically  stable 
at  octahedral  geometries. 

We  believe  that  electronic  stability  could  be  further  increased 
by  arranging  a  closed-shell,  ionicity-enhancing  structure  for  multiply 
charged  anions.  This  can  be  achieved  by  a  proper  choice  of  the 
central  atom.  Indeed,  recent  theoretical  studies  on  thermochemical 
properties  of  doubly-charged  hexahalogenometallates  involving  Ti, 
Zr,  and  Hf  demonstrated  their  significant  stability  with  respect  to  the 
metal  atom  plus  three  halogen  molecules.^^  These  promising  results 
were  obtained,  however,  at  low  theoretical  levels,  motivating  us  to 
extend  the  investigations  to:  a)  higher  levels  of  theory,  b)  analogous 
closed-shell  singly-  and  triply-charged  systems. 

We  choose  highly  electropositive,  bulky  transition  metals  such 
as  Ta,  Zr,  and  La  for  the  singly-,  doubly-,  and  triply-charged  anions, 
respectively,  to  ensure  highly  ionic  character  in  the  coordinated 
complex  and  to  make  the  distance  between  the  negatively  charged 
ligands  as  large  as  possible.  Hf  would  be  more  consistent  with  the 
choice  of  Ta  and  La,  but  it  is  well  known  that  chemical  properties  and 
bond  distances  are  practically  the  same  for  the  hafnium  and 
zirconium  atoms.3 1 

Synunetry  properties  and  orbital  energies  of  the  metal  M  and 
ligand  L  allow  us  to  make  an  educated  guess  about  the  properties  of 
the  molecular  orbitals  in  MLa**'.  In  Table  1  we  show  the  symmetry- 
adapted  orbitals  (SAO)  resulting  from  the  valence  atomic  orbitals 
(AO),  as  well  as  orbital  energies  in  M  and  L.  The  orbital  energies 
permit  prediction  of  the  ordering  and  nature  of  the  molecular 
orbitals  in  the  ML6°'  complexes.  In  particular,  it  is  clear  that  high- 
energy,  occupied  molecular  orbitals  will  be  dominated  by  the  ligands' 
np-typc  AO’s.  Mixing  with  the  M  orbitals  will  be  largest  for  Ta  and 
smallest  for  La.  The  ligands'  np  orbitals  may  be  separated  into  two 
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groups:  parallel  (np||)  and  perpendicular  (npl)  to  the  ML  bond.  The 
npll  group  contains  orbitals  of  the  same  symmetry  as  the  valence 
orbitals  of  the  M  center,  hence  bonding  interaction  can  develop.  On 
the  other  hand,  in  the  npl  group  only  the  tiu  orbital  can  interact, 
presumably  in  an  antibonding  manner,  with  the  M  valence 
counterpart.  All  of  the  other  ligand  orbitals  are  nonbonding.  We 
conclude  that  a  few  highest  occupied  orbitals  of  the  closed-shell 
MXgn-complex  will  possess  either  a  nonbonding  (tig,  t2u.  t2g)  or 
antibonding  (tiu)  character  and  that  the  doublet  ground  state  of  the 
system  will  be  triply  degenerate  at  the  Ojj  geometry. 

Based  on  such  an  analysis,  we  selected  for  our  calculations  the 
following  triply  charged  anions:  A1F63*,  GaF63*,  InFe^*,  TIF6^‘,  LaF63- 
and  LaCl6^*  as  well  as  the  doubly  charged  ZrFg^-  and  singly  charged 
TaFe*.  All  of  these  systems  satisfy  the  MFk+s^'.  MFk+2^'  and  MFk+r 
rules  (see  above). 

III.  Computational  Aspects 

A.  Transition  Metal  Anions. 

For  LaF63*.  LaCl63-,  ZrFa^*,  TaFa",  LaFsZ-,  LaCls^-,  ZrFs*,  TaFs. 
LaFa",  LaCU*,  ZrFa  and  TaF4+  species,  we  used  the  relativistic,  compact 
effective  potentials  of  Stevens  at  al.  together  with  their  valence, 
energy-optimized,  shared-exponent,  contracted  gaussian  basis  sets, 
which  we  label  SBKJ.32  For  the  M  atoms,  the  (n-l)s,  (n-l)p,  ns,  np, 
and  (n-l)d  orbitals  are  chemically  active  and  take  advantage  of  the 
double-zeta  sp  and  triple  zeta  d  basis  sets.  For  the  L  atoms,  only  the 
valence  ns  and  np  orbitals  are  chemically  active,  and  they  are 
described  within  a  double-zeta  quality  basis  set. 

Next,  the  SBKJ  basis  sets  were  supplemented  with  different  sets 
of  diffuse  functions.  The  resulting  differences  in  molecular 
geometries  and  Koopmans’  theorem  predictions  of  electron 
detachment  energies  were  carefully  monitored.  For  each  M  we  tested 
one-  and  for  each  L  one-  and  two-term  sets  of  diffuse  functions.  The 
exponents  of  the  one-term  diffuse  basis  sets  were  taken  to  be  equal 
to  1/3  of  the  lowest  exponent  in  the  original  basis.  For  the  two-term 


diffuse  sets,  a  geometrical  progression  was  assumed  with  the 
progression  constant  equal  to  2.4.  Numerical  tests  described  in  the 
next  Section  suggested  use  of  the  original  SBKJ  basis  set 
supplemented  with  one  set  of  diffuse  functions  on  the  L's  only.  The 
resulting  basis  set  for  ML6"'  complexes,  labeled  SBKJ+diff,  consists  of 
106  contracted  gaussian  basis  functions.  It  should  be  stressed  that 
due  to  the  absence  of  polarization  functions,  our  electron  detachment 
energies  are  expected  to  be  underestimated. 

We  performed  geometry  optimization  for  the  octahedral  ML6"' 
complexes  with  the  SBKJ  and  SBKJ+diff  basis  sets  using  the  GAMESS 
program  package. ^3  Gradient  geometry  optimization^^  was  followed 
by  calculation  of  the  numerical  hessian  matrix.  Normal  mode 
analyses  were  performed  and  IR  intensities  were  then  calculated. 

We  calculated  two  quantities  that  characterize  the  charge 
distribution  in  the  multiply  charged  anions.  First,  the  radial  extent  of 
the  distribution  is  characterized  by  the  mean  value  of  (<R2>). 
Second,  the  effective  charges  assigned  to  the  M  and  L  atoms  were 
determined  from  the  fit  to  the  electrostatic  potential  at  points 
selected  according  to  the  Merz-Singh-Kollman  scheme.^^  For  atomic 
radii,  we  used  1.34,  1.45,  1.69  A  for  Ta,  Zr,  and  La,  respectively ;3 1  for 
the  halogens  we  used  the  default  values  (1.35  A  (F)  and  1.70  A  (Cl)) 
.3  6  These  calculations  were  performed  using  the  Gaussian  92 
package36  with  the  SBKJ+diff.  basis. 

Vertical  detachment  energies  were  calculated  directly,  using 
the  electron  propagator  theory  (EPT)  approach  [37-42],  and 
indirectly,  i.e.,  calculating  the  energy  of  the  n-electron  closed-shell 
anion  and  (n-l)-electron  open-shell  species  using  commonly 
available  Hartree-Fock  (HF)  and  the  Mpller-Plesset  perturbation 
theory  (MP)  mcthods.33.36 

In  the  EPT  approach,  two  types  of  self-energy  approximations 
are  employed.  In  the  first,  second  order  poles  are  calculated;  these 
results,  designated  by  the  abbreviation  EPT2,  typically  exaggerate 
electron  correlation  effects.  The  second  method  includes  all  third 
order  terms  and  many  infinite  order  terms  as  well.  Energy- 
independent  terms  through  third  order  are  added  to  an  energy- 
dependent  self-energy  that  includes  all  ring  and  ladder  diagrams 


through  all  orders  in  the  fluctuation  potential.^^  This  method  is 
therefore  identified  by  the  abbreviation  EPT3+.  Third-order  methods 
have  been  thoroughly  examined  for  electron  affinities.‘*0.44  while  the 
fluoride  anion  requires  additional  terms  in  order  to  obtain  an 
accurate  detachment  energy,^'*  third  order  methods  are  generally 
reliable  for  fluorine-containing  compounds. The  usual  pole-search 
algorithms  for  second  order  have  been  used. 39  For  EPT3+  calculations, 
a  ne>v  algorithm,  based  on  eigenvalue  procedures  used  in 
configuration  interaction  calculations,  has  been  employed."^^  A 
modified  version  of  the  EPT90  code^6  was  used  to  perform  these 
calculations. 

Indirect  calculations  of  VDE's  were  restricted  to  the  electronic 
state  of  obtained  by  the  removal  of  an  electron  from  the 

highest  occupied  molecular  orbital  (HOMO)  of  ML6°*.  The  results 
discussed  in  the  next  Section  demonstrate  that  the  ordering  of  the 
electronic  states  of  predicted  on  the  basis  of  the  Koopmans' 

theorem  remains  the  same  in  the  EPT  calculations. 

We  started  from  one-configuration  calculations  to  determine 
the  relaxation  correction  to  the  Koopmans’  theorem  prediction. 
Calculations  for  the  triply -degenerate  doublet  state  of  the  neutral 
were  performed  in  two  ways.  First,  a  symmetry-unbroken  solution 
was  obtained  with  the  hole  delocalized  over  six  ligands.  This  was 
achieved  in  the  state-averaged  (SA)  SCF  calculation^'^  with  three 
configuration  state  functions  involved.  It  has  been  reported, 
however,  that  single  configuration  calculations  with  a  core'^^  or 
valence^^  hole  delocalized  over  equivalent  atoms  may  produce 
unphysical  electron  binding  energies.  More  reliable  results  were 
obtained  when  the  hole  was  localized  in  one  configuration  calculation, 
i.e.,  when  a  symmetry  broken  solution  was  produced.'***'*^  Hence,  we 
reinforced  the  SA  SCF  calculation  with  UHF  calculations  in  which  we 
allowed  for  localization  of  the  hole  on  on*  ligand.  Of  course,  after 
inclusion  of  correlation  effects,  both  single-configuration  approaches 
should  lead  to  the  same  final  result.  Our  correlated  MPn  (n=:2,3,4) 
calculations  were  performed  with  the  UHF  reference  function.3  6 
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B.  Nontransion  Element  Hexafluorides 

Bond  lengths  and  harmonic  frequencies  for  the  isoelectronic 
AIF63‘,  GaF63-,  lnF6^'  and  T1F6^'  as  well  as  for  AIF3,  GaFs,  InFs,  TIF3. 
A1F4‘,  GaFa*,  InFa-  and  T1F4-  have  been  obtained  using  RHF 
approximation  and  the  Los  Alamos  Pseudopotentials, 50  which  we 
label  LAP.  We  used  Dunning's  valence  double-zeta  basis  set 
(9s5p/3s2p)  for  the  fluorine  atoms^i  and  a  valence  2s2p  basis  sets 
for  the  aluminum,  gallium,  indium  and  tallium,50  extended  by 
polarization  d-functions  on  all  atoms  (LAP/DZ+d)  and  diffuse  s-  and 
p-functions  on  the  fluorine  atoms  (LAP/DZ+dfdiff.).  Exponents  for  the 
polarization  d-functions  0.198  (Al),  0.207  (Ga),  0.160  (In)  and  0.146 
(Tl)  were  taken  from  Ref.  52  and  exponents  of  the  diffuse  s-  and  p- 
functions  0.074  were  taken  from  Ref.  53.  All  calculations  with  Los 
Alamos  pseudopotential  were  performed  using  the  Gaussian  92 

package. 36 


IV.  Results 

A.  Nontransition  metai  hexafluorides. 

The  results  of  SCF  (at  LAP/DZ+d  and  LAP/DZ+d+diff.)  geometry 
optimization  and  frequency  calculations  are  presented  for  AIF3,  GaFs, 
InF3  and  TIF3  in  Table  IIA,  for  A1F4',  GaF4",  InF4'  and  TIF4*  in 
Table  IIB  and  for  AlFe^*,  GaFe^-,  inFe^*  and  TlFe^-  in  Table  IIC. 
Experimental  data  are  available  for  AIF3,  GaF3,  AIF4'  and  for  all 
triply  charged  hexafluorides.  Our  optimized  bond  lengths  and 
frequencies  are  nearly  the  same  for  the  LAP/DZ+d  and 
LAP/DZ+d+diff.  calculations,  and  agree  with  available  experimental 
data  within  0.05  A  for  bond  lengths  and  30  cm*l  for  vibrational 
frequencies  (Tables  IIA  and  IIB). 

The  optimized  bond  lengths  R(A1-F)=:  1.89- 1.91  A  of  A1F6^’, 
R(Ga-F)=  1.92- 1.96  A  of  GaF63-,  R(In-F)=2.05-2.06  A  of  InFe^-  and 
R(Tl-F)=2.22-2.24  A  of  TlF63-can  be  compared  with  the  experimental 
data  for  these  anions  in  crystal  salts  R(A1-F)=1.81  A  in  Na3AlF6,54 


Cs2NaAlF6;55  R(Ga-F)=1.89  A  in  (NH4)3GaF6;56  R(In-F)=2.04  A  in 
(NH4)3lnF657  and  R(Tl-F)=l.96  A  in  Rb3TlF6.57  From  these  data  one 
can  see  that  the  difference  between  the  calculated  and  corresponding 
experimental  bond  lengths  is  high  for  A1F6^*  (0.1  A),  decreases  in 
GaF63'  (0,03-0.07  A)  and  further  decreases  for  InF6^*  (0.01 -0.02A), 
but  the  difference  is  very  large  in  T1F6^'  (0.26-0.28  A).  Why  our 
bond  length  is  too  high  for  TlFg^-  js  ^^t  clear  at  this  time.  The 
experimental  bond  lengths  increase  from  AIF3  to  A1F6^'  and  from 
GaF3  to  GaF6^'  by  approximately  0.18  A,  and  the  same  variation  is 
found  in  our  calculations.  From  this  point  of  view,  the  calculated 
difference  in  bond  lengths  (0.24  A)  for  TIF3  and  T1F6^’  may  not  be 
unreasonable. 

All  of  AlFe^',  InFg^-  and  T1F63-  have  local  minima  at  the 

octahedral  structure,  and  thus  all  of  the  computed  vibrational 
frequencies  are  real  (see  Table  IIC).  Therefore,  if  we  keep  the 
number  of  electrons  constant,  these  trianions  are  geometrically 
stable.  However,  the  disagreement  between  our  calculated  and  the 
corresponding  vibrational  frequencies  of  crystal  salts  is  unexpectedly 
high  (80-160  cm'l).  At  the  same  level  of  theory  the  maximum 
difference  between  experimental  and  calculated  frequencies  for 
AIF3,  A1F4‘  and  GaF3  is  only  30  cm- 1  (see  Tables  IIA  and  IIB). 
Moreover,  when  we  improved  the  basis  sets  and  included  diffuse 
functions,  the  disagreement  became  even  worse. 

We  suspect  that  these  large  differences  may  be  related  to  the 
electronic  instability  of  the  anions.  In  particularly,  the  MXa  in  the 
solid  state  would  probably  bear  a  smaller  excessive  charge  than  that 
assumed  in  our  ab  initio  calculations.  The  smaller  excess  charge 
would  make  the  species  less  floppy,  in  agreement  with  the 
experimental/theoretical  difference.  Moreover,  the  electronic 
instability  may  enhance  the  coupling  between  the  vibrations  of  MX6 
and  the  neighboring  species.  Therefore,  some  vibrations  should  be 
assigned  to  the  isolated  anions,  but  rather  to  a  collective  motion  also 
involving  the  cations  and  even  the  lattice  as  a  whole.  If  so,  the  so 
called  "experimental  force  fields"58.59  for  these  anions  should  be 
revised. 


The  orbital  energies  of  the  tig-HOMO  are  positive  for  A1F63- 
and  GaFg^-  and  therefore  both  these  anions  are  certainly 
electronically  unstable,  because,  according  to  our  experience,  when 
electron  relaxation  and  electron  correlation  are  taken  into  account 
the  electronic  stability  of  multiply  charged  anions  decreases.  InFe^- 
and  T1F63-  both  have  negative  orbital  energies  of  the  tig- HO  MO: 
-0.41  eV  and  -0.20  eV,  respectively.  Therefore,  both  the  anions  are 
electronically  stable  at  the  Koopmans  level  of  theory.  However, 
previous  ab  initio  calculation  for  BeF42*,12  MX32-,  where  M=Li-K  and 
X=F-C1, 13.14  Se042-  and  Te042-.l8  have  shown  that  electron  correlation 
and  electron  relaxation  corrections  decrease  electronic  stability  by  1- 
2  eV.  Therefore  we  conclude  that  all  free  AlFe^*,  GaF63",  InF63-  and 
TlFe^-  are  not  electronically  stable  species. 

B.  Transition  metal  hexafluorides. 

The  results  of  SCF  geometry  optimization  for  TaF6', 

LaF63‘,  and  LaCl6^'  closed-shell  anions  with  the  SBKJ  basis  set  are 
presented  in  Table  HD.  We  also  report  vibrational  frequencies  (all 
real)  and  IR  intensities  for  these  species.  We  have  demonstrated  that 
augmenting  the  SBKJ  set  with  diffuse  functions  on  X  and/or  M  affects 
the  M-L  equilibrium  distance  by  less  than  0.01  A.  Results  of 
Cederbaum  and  collaborators  indicate  that  correlation  effects  barely 
affect  equilibrium  geometries  of  tetrahedral  and  triangular  doubly 
charged  anions. Hence,  our  further  calculations  were  performed 
at  the  SBKJ  equilibrium  geometries. 

Calculated  geometries  and  vibrational  frequencies  for  isolated 
TaF6'  and  ZrF6^'  anion  are  in  good  agreement  with  available 
experimental  data  (see  Table  IID)  from  the  solid  phase.  The 
discrepancies  in  frequencies  are  less  than  25  cm*l  for  TaF6*.  For 
ZrF62*,  the  SCF/SBKJ  frequencies  are  systematically  lower  than  the 
solid  state  data.^O-^*  but  the  discrepancies  do  not  exceed  45  cm-l.  For 
triply  charged  LaFe^-  and  LaCl6^"  species,  however,  discrepancies 
are  more  significant  (  108  cm"l  for  valent  symmetrical  vibration. 
Table  IID).  These  triply  charged  anions  are  electronically  barely 
stable  or  even  unstable  (see  below).  This  may  cause  larger 


disagreement  with  experimental  data  (see  above  discussion  for 
nontransition  metal  hexafluorides). 

The  values  of  <R2>  presented  in  Table  II  demonstrate  how  the 
spatial  extent  of  the  electronic  charge  distribution  changes  in  the 
sequence  of  multiply  charged  but  otherwise  similar  anions.  For 
comparison,  the  values  of  <R2>  for  F*  and  Cl*  are  15.8  and  36.9  a.u., 
respectively.  As  expected,  the  spatial  extent  is  particularly  large  in 
the  case  of  three  excess  electrons.  Interestingly,  inspection  of  the 
orbitals  in  the  HOMO  region  indicates  that  they  are  not  dominated  by 
the  most  diffuse  functions  which  are  available  in  the  basis  set.  We 
interpret  this  feature  as  the  first  indication  that  the  triply  charged 
anion  may  be  electronically  stable.  Effective  atomic  charge,  presented 
in  the  fourth  column  of  Table  II,  reflects  the  increasing  ionic 
character  in  the  Ta,  Zr,  La  sequence.  Replacing  Fs  by  Cl's  increases 
substantially  the  spatial  extent  of  the  electronic  charge  distribution 
and  slightly  decreases  the  ionicity  of  the  complex.  The  reported 
values  of  effective  atomic  charges  allow  consideration  of  ZrF62-, 
LaFe^*,  and  LaCle^*  as  ionic  complexes  with  six  singly  charged  ligands 
coordinated  to  the  M(6-n)+  core,  consistent  with  the  ionic  model  of 
Cederbaum  and  collaborators. *3. 14 

Vertical  detachment  energies  from  the  HOMO  orbital  of  the 
closed-shell  species  are  reported  in  Table  III.  All  the  species  are 
electronically  stable  at  the  Koopmans  theorem  level.  Augmenting  the 
SBKJ  basis  set  with  one  set  of  L(sp)  diffuse  functions  centered  on  L's 
may  change  the  HOMO  orbital  energy  even  by  0.4  eV  (LaFe^'  case). 
However,  further  extensions  of  the  basis  set  with  two  sets  of  ligand 
diffuse  functions  or  one  set  of  M  diffuse  functions  changes  the  HOMO 
orbital  energy  by  at  most  a  few  hundredths  of  an  eV.  This  further 
supports  our  choice  of  SBKJ+diff  as  the  basis  set  for  the  production 
runs. 

The  SCF  relaxation  energies  for  F*  and  Cl*  calculated  within  the 
SA  SCF  and  UHF  procedures  are  very  large  but  quite  similar.  For 
MF6<*',  however,  the  SA  SCF  relaxation  energies  do  not  exceed  0.65 
eV.  At  the  same  time  the  UHF  relaxation  energies  may  be  as  large  as 
3.51  eV  and  the  triply  charged  anions  are  predicted  to  be  unstable  at 
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this  level  of  theory.  Interestingly,  the  relaxation  effects  are  much 
smaller  for  LaCl6^',  as  they  are  for  Cl-  compared  to  F*. 

There  is  a  significant  difference  in  the  magnitude  of  the 
relaxation  contribution  to  the  VDE  when  the  UHF  or  SA  SCF  scheme 
are  applied,  see  also  Refs.  48,  49  and  14b.  In  all  of  the  fluoride- 
containing  species  studied  here,  the  UHF  based  VDE  is  further  from 
the  correlated  result  than  the  SA  SCF  prediction;  the  later  however, 
systematically  overestimates  the  magnitude  of  the  binding  energy  in 
all  MLe^-  systems.  Hence,  neither  the  symmetry-broken  UHF  nor  the 
SA  SCF  delivers  reliable  VDE’s. 

Inclusion  of  correlation  effects  at  the  Mdller-Plesset  level  of 
theory  increases  the  detachment  energy  in  comparison  with  the  UHF 
predictions.  The  convergence  in  consecutive  orders  of  the 
perturbation  theory  is,  however,  not  uniform,  with  the  MP2  and  MP4 
results  being  usually  close  to  each  other.  The  similarity  between  the 
MP2  and  MP4  VDE's  should  not  be  considered  as  a  firm  argument 
that  the  converged  result  has  been  reached.  For  instance,  for  CN*,  the 
MP2  and  MP4  results  are  quite  similar  but  inaccurate.^ 5  The 
experimental  electron  detachment  energies  for  F*  and  Cl*  arc  3.40 
and  3.62  eV,  respectively .^2  Therefore,  our  MP4  results  for  these 
species  are  off  by  0.28  and  0.50  eV,  respectively.  Hence,  the 
electronic  instability  of  LaF6^'  and  LaCl6^*  predicted  at  the  MP4  level 
is  probably  misleading. 

Vertical  detachment  energies  calculated  within  the  EPT  method 
also  exhibit  a  non-uniform  convergence.  The  EPT2  correction  to  the 
Koopmans  theorem  result  is  destabilizing  and  overshoots.  Triply- 
charged  anions  are  unstable  at  this  level  of  theory.  The  EPT3+  results 
are,  however,  much  more  positive,  predicting  LaFe^-  to  be 
electronically  stable  and  LaCl63’  marginally  unstable.  Since  the  EPT3+ 
result  for  Cl*  is  underestimated  by  0.21  eV,  we  think  that  LaCle^" 
may  be  electronically  stable.  Contrariwise,  the  EPT3+  result  for  F* 
puts  the  electronic  stability  of  LaF63-  in  question.  The  slow 
convergence  of  EPT  methods  in  the  case  of  F*  is  well  documented.^'* 
An  encouraging  observation  is  that  EPT  VDE's  of  the  fluorine - 
containing  compounds  are  usually  reliable.^^ 
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Our  MP4  and  EPT3+  results  demonstrate  that  TaF6'  is  a  striking 
example  of  a  superhalogen,  i.e.,  an  anion  with  the  electron  binding 
energy  exceeding  that  of  Cl*  ion.2l.63.64  its  EPT3+  value  of  VDE  (10.61 
eV)  is  at  least  as  large  as  the  experimental  result  of  10  ±  0.5  eV  for 
the  "unusual"  AuFa’  species.27,28  Previously  calculated  VDE(TaF6* 
)-8.4  eV  by  DV-Xa  method20a  is  somewhat  underestimated. 

Interestingly,  in  the  doubly-charged  molecular  anion  ZrF62-,  the 
second  excess  electron  is  bound  more  strongly  than  the  first  excess 
electron  in  the  leading  monoatomic  anion  Cl'.  This  conclusion  is  valid 
both  at  the  MP4  and  EPT3+  level.  The  value  of  4.99  eV  (EPT3+  result) 
is  comparable  with  the  estimation  (4.9-5. 1  eV)  of  VDE  in  TeF82-, 
claimed  to  be  the  MXn2'  dianion  with  the  largest  VDE.i5  Since  many 
ionic  solids  containing  ZrF62-  dianion  are  known,39a  one  could  try  to 
produce  this  unusual  species  in  sputtering  experiments. 

As  expected  from  the  symmetry  analysis  and  orbital  energies 
presented  in  Table  1,  for  all  MX6‘^~  species  the  three  highest  occupied 
molecular  orbitals  are  of  tiu,  tig,  and  t2u  symmetry.  In  Table  IV, 
VDE's  for  the  transitions  to  the  doublet  Tiu»  Tig,  and  T2u  states  are 
reported  for  all  fluoride-containing  complexes.  The  ordering  of  these 
states  depends  on  the  system,  but  for  every  species  it  remains  the 
same  at  the  KT  and  EPT3+  level.  More  interestingly,  the  separation 
between  the  first  and  third  state  predicted  at  the  KT  level  remains 
within  0.04  eV  of  the  separation  at  the  EPT3+  level.  Moreover,  the 
pole  strengths  in  the  EPTS-t-  calculations,  which  are  reported  in  Table 
IV,  never  drop  below  0.91.  These  two  factors  support  our  claim  that 
not  only  doubly-  but  also  triply-charged  octahedral  anions  may  be 
electronically  stable. 

Finally,  we  briefly  discuss  thermodynamic  stability  of  the 
ML 6*^'  species  with  respect  to  two  dissociation  steps: 

ML6“-  -->  ML5(“-*)-  +  L'  (la) 

ML6"-  ->  ML4(“-2)-  +  2L-  (lb) 

Equilibrium  geometries,  vibrational  frequencies,  and  IR 
intensities  for  the  product  molecules  ML5(“*l)*  (Dsh)  and  ML4("'2)- 


(Td)  were  determined  at  the  SCF/SBKJ  level  and  are  reported  in  Table 
V. 

Energy  changes  for  the  (la)  and  (lb)  reactions  were 
determined  at  the  MP2  level  with  the  SBKJ+diff  basis  sets  and  are 
reported  in  Table  VI.  In  these  MP2  calculations,  the  four  lowest 
occupied  orbitals  were  not  correlated.  The  results  demonstrate  that 
the  decomposition  of  ML63‘  will  be  exothermic  in  both  reaction  steps. 
The  amount  of  energy  released,  533  and  475  kJ/mol  for  LaFe^-  and 
LaCl6^%  respectively,  is  significant  and  the  entropy  factor  Id 
increase  exothermicity.  Our  results  also  indicate  that  LaFs  id 
LaCls^'  species  are  promising  candidates  for  electronically  and 
geometrically  stable  penta-coordinate  dianions.  Their  KT  values  of 
VDE  are  5.4  and  4.2  eV,  respectively. 

The  la  decomposition  of  ZtFs^'  is  predicted  to  be  slightly 
exothermic.  The  product  ZrFs*  is  thermodynamically  stable.  In 
addition,  ZrFy  displays  properties  of  a  remarkable  penta-coordinate 
superhalogen.  Its  VDE  at  the  KT  level  is  11.4  eV. 

Finally,  TaF6%  which  was  recognized  above  as  a  remarkable 
hexa-coordinated  superhalogen,  is  thermodynamically  stable  with 
respect  to  the  decompositions  la  and  lb. 

V.  Summary 

We  have  studied  local  geometric  and  electronic  stability  of  the 
closed-shell  octahedral  species  AIFe^-,  GaFg^-,  InFa^-,  TlFe^-,  TaFe*, 
ZrFe^",  LaFe^-,  and  LaCle^'  using  ab  initio  methods. 

We  have  found  that  local  minima  exist  for  these  species  at 
octahedral  geometries. 

AlFe^'  and  GaFe^*  ions  have  positive  HOMO  energies  and 
therefore  are  not  electronically  stable  even  at  the  Koopmans  theorem 
level.  InFe^*  and  TlFe^*  have  negative  orbital  energies  of  -0.41  eV 
and  -0.21  eV.  Because  both  energies  are  low  ^nd  correction  from 
electron  correlation  and  electron  relaxation  should  destabilize  both 
anions  by  1-2  eV,  we  concluded  that  all  four  of  these  triply  charged 
anions  are  not  electronically  stable  species. 


Calculated  vibration  frequencies  for  electronically  unstable 
A1F63-,  GaF63-,  InF63-  and  T1F63-  species  are  underestimated  by 
100-160  cm‘^,  while  for  electronically  stable  species  (AIF3,  A1F4-, 
GaF3,  ZrF4  and  ZrF6^'),  these  numbers  differ  from  experimental  data 
by  no  more  than  45  cm'  ^ .  We  consider  these  larger  disagreements 
for  triply  charged  anions  to  be  a  result  of  a)  different  excess  negative 
charges  operative  in  the  solid  state  and  b)  essential  interactions 
between  vibrational  modes  of  isolated  anions  and  the  lattice  (cations) 
for  electronically  unstable  anions.  Therefore,  experimental 
vibrational  spectra  for  solids  containing  multiply  charged 
electronically  unstable  anions  should  be  treated  by  taking  collective 
modes,  including  cations  and  anions,  into  account. 

Electronic  stability  of  TaF6*,  LaF6^',  and  LaCl6^'  was 

investigated  using  the  electron  propagator  and  the  M0ller-Plesset 
perturbation  approach.  Our  results  indicate  that  triply-charged 
molecular  anions  may  be  electronically  stable  as  isolated  species.  The 
convergence  of  the  EPT  and  MP  series  was  found  to  be  slow  but  the 
error  analysis  for  isolated  F*  and  Cl*  supports  the  hypothesis  of 
electronic  stability  of  ML^^'  species  studied  here. 

Both  triply-charged  species  were  found  to  be 
thermodynamically  unstable  with  respect  to  two  consecutive 
unimolecular  decompositions.  The  electronic  energy  released  would 
be  533  and  475  kJ/mol  for  LaF^^'  and  LaCl6^',  respectively.  Hence 
they  are  promising  components  of  the  high  energy  density 
propellants. 

More  advanced  ab  initio  calculations  as  well  as  exploration  of 
geometries  probed  in  the  zero-point  vibrational  motion  of  ML63* 
species  are  required  to  make  conclusions  about  the  long-lived 
metastability  of  these  interesting  species. 

The  ZrF6^'  species  was  found  to  be  a  promising  candidate  for  an 
extremely  stable  molecular  dianion.  Its  VDE  value  amounts  to  4.99 
eV  at  the  EPT3+  level  of  theory  and  relatively  low  thermodynamic 
instability  with  respect  to  ZrFs*  +  F*  supports  our  claim. 

The  TaFfi"  species  was  recognized  as  a  thermodynamically 
stable  superhalogen  with  the  VDE  value  of  10.6  eV  at  the  EPT3+  level 
of  theory. 
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Table  I.  Symmetry-adapted  orbitals  (SAO)  that  can  be  formed  from 
the  ligand  (L)  and  metal  (M)  atomic  orbitals  (AO)  for  octahedral 
complexes  ML6"‘.  Orbital  energies  of  L  and  M  in  a.u.  are  calculated 
with  the  SBKJ+diff  basis  set.  Symbols  np||  and  npl  denote  ligands'  np 
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Table  IIB.  SCF  optimized  central  atom-ligand  distances  (Rml  in  A),  vibrational  frequencies  (cnrO. 
and  IR  intensities  (D^/(A2  amu);  in  parenthesis)  for  the  octahedral  closed-shell  anions.  Spatial 
extent  of  the  SCF  electronic  charge  distribution  <R2>  in  a.u.  Effective  atomic  charge  (Mulliken)  on 
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Table  IlC.  SCF  optimized  central  atom-ligand  distances  (Rml  in  A),  vibrational  frequencies  (cm  *), 
and  IR  intensities  (D^/CA^  amu);  in  parenthesis)  for  the  octahedral  closed-shell  anions.  Spatial 
extent  of  the  SCF  electronic  charge  distribution  <R^>  in  a.u.  Effective  atomic  charge  (Mulliken)  on 
the  ligand  atom  (Ql)  in  a.u. 
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Table  IID.  SCF  optimized  central  atom-ligand  distances  (Rml  A),  vibrational  traiiicncics  (cm 
and  IR  intensities  (D2/(A2  amu);  in  parenthesis)  tor  the  octahedral  closed  shell  amons.  Spatial 
extent  of  the  SCF  electronic  charge  distribution  <R2>  in  a.u.  Effective  atomic  charge  on  the  ligand 
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Table  III.  Electron  vertical  detachment  energies  (in  eV)  from  the 
highest  occupied  molecular  orbital  of  the  closed-shell  anion. 
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Final  state 
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Tabic  Va.  SCI-VSBKJ  oplimized  metal -ligand  axial  (Rax)  and  equatorial  (Req)  distances  (in  A), 
vibrational  frequencies  (cm  '),  and  IR  intensities  (D2/(A2  amu);  in  parenthesis)  for  the  D^j,  closed 
shell  complexes  MLs- 
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Data  from  Ref. 


Table  IV,  EPT  electron  vertical  detachment  energies  (in  eV)  from  the 
closed-shell  *Aig  octahedral  anion  to  different  doublet  states.  Values 
of  the  EPT3'*’  pole  strengths  in  parenthesis. 


Table  Vb.  SCF/SDKJ  optimized  metal-ligand  distance  (in  A),  vibrational  frequencies  (cm'l),  and  IR 
intensities  (D2/(A2  amu);  in  parenthesis)  for  the  Td  closed-shell  complexes  ML4 
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Table  VI.  MP2  electronic  energy  changes  in  kJ/mol  for  the  reactions 
la  and  lb.  The  zero  point  vibrational  corrections  are  not  included. 
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